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ABSTRACT 


One -dimensional,  unsteady  flow  theory  of  an  ideal 
gas  is  reviewed  as  it  applies  to  guns  and  ballistic  model 
launchers.  The  method  of  characteristics  is  developed 
for  the  general  case  including  chambrage  and  finite 
chamber  volume.  Results  for  helium  as  a  propellant 
are  given  from  which  the  required  chamber  size  can  be 
determined.  The  launch  velocity  for  combinations  of 
launch  tube  and  chamber  geometries  is  computed  assum¬ 
ing  no  friction  nor  heat  loss  and  an  evacuated  bore. 
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NOMENCLATURE 


A  Cross-section  area 

a  Acoustic  speed 

m  Mass 

p  Pressure 

S  Entropy 

8  Distance 

T  Temperature 

t  Time 

u  Velocity 

uo  Reference  launch  velocity  (Ai/A,  =  1,  xch  ^ 

Ve  Volume  of  chamber 

Vl,  Volume  of  launch  tube 

xch  Chamber  length 

y  Ratio  of  specific  heats 

p  ■■  Density 


SUBSCRIPTS 


c  Chamber  initial  value 

cli  Conditions  at  end  of  chamber 

F  Forward  face  of  piston 

L  Sonic  exit  limit;  launch  tube 

o  Standard 

p  Piston  or  projectile 

R  Rear  face  of  piston 

t,  s  Partial  derivatives 

1  Upstream  of  chamber  exit  |  chamber 

2  Downstream  of  chamber  exit  J 


7 


AEDC-TR-61-1 


t 

P=  Al 

“M  «c 

t 

11 

uAc 

a 

a/ac 

P 

P  C«  ®c 

P/Pc 

Initial  chamber  conditions 

mM 

Projectile  mass 

8 


AEDC-TR-61-1 


INTRODUCTION 


The  light-gas  gun  type  of  ballistic  model  launcher  has  made  it  pos¬ 
sible  to  simulate  flight  environment  at  very  high  speeds  more  closely 
than  by  any  other  ground  experimental  facility.  A  light  gas  (hydrogen 
or  helium)  is  necessary  as  a  propellant  because  the  associated  high 
acoustic  speed  allows  the  chamber  pressure  to  be  transmitted  to  the 
projectile  base  with  lesser  decrease  at  high  speeds  than  when  a  heavier 
propellant  is  used. 

The  highest  velocities  have  been  achieved  with  launchers  in  which 
the  propellant  light  gas  is  raised  to  a  final  high  energy  state  by  com¬ 
pression  with  a  free  piston.  Velocities  in  the  range  between  25,000  and 
30,  000  ft/sec  are  regularly  obtained  using  light  projectiles. 

One  of  the  important  aspects  of  the  design  of  launchers  is  the  effect 
of  geometrical  variables  of  the  configuration  on  the  performance.  The 
launch  tube  length,  chamber  length,  and  chambrage  (chamber  cross- 
section  to  launch  tube  cross-section  area  ratio)  are  the  essential  dimen¬ 
sions.  The  calculation  of  the  unsteady  flow  in  one  dimension  by  the 
method  of  characteristics  appeared  to  be  the  most  practicable  approach 
to  the  problem.  Reference  1  provides  a  general  background  for  the 
analysis,  and  Seigal  (Ref.  2)  treats  the  effect  of  chambrage  for  infinitely 
long  chambers.  The  effect  of  finite  chamber  length  is  considered  by 
Heybey  (Ref.  3)  for  a  launch  tube  the  same  diameter  as  the  chamber. 
Calculations  for  ratios  of  specific  heats,  y ,  of  1.  25  and  1.  4  are  shown 
in  Ref.  2,  and  the  method  of  calculation  is  given  in  Ref.  3,  It  was  con¬ 
sidered  most  desirable  to  make  the  characteristics  calculations  for  fmite 
chamber  volume,  chambrage,  and  y  -  1.66  corresponding  to  helium  in 
order  to  have  useful  results  for  light -gas  launchers  that  can  be  verified 
experimentally.  In  a  two -stage  configuration,  the  volume  after  compres¬ 
sion  is  likely  to  become  very  small,  and  it  varies  widely  with  operating 
conditions.  It  is  therefore  obvious  that  the  compression  tube  must  be 
properly  matched  to  the  launch  tube  to  obtain  the  optimum  performance 
of  the  system.  The  general  characteristics  solutions  provide  the  infor¬ 
mation  required  for  finding  the  minimum  chamber  for  a  fixed  launch  tube. 

The  characteristics  results  provide  a  starting  point  for  estimating 
the  effect  of  heat  losses  from  the  propellant.  Since  the  method  la  a  step- 
by-step  process,  it  lends  itself  readily  to  the  introduction  of  projectile 
friction,  real  gas  properties,  and  relaxation  times . 


Manuscript  released  by  author  February  1961. 
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UNSTEADY,  ISENTROPIC  PROCESSES  IN  ONE-DIMENSIONAL  FLOW 


GENERAL  RELATIONS 


The  concept  of  one -dimensional  flow  treated  here  corresponds  to  the 
motion  of  a  gas  in  a  constant  area  tube  in  which  the  properties  are  uniform 
across  any  section.  No  heat  is  transferred  across  the  tube  walls,  and 
viscosity  is  neglected. 

Basic  equations  for  one -dimensional  flow  of  an  ideal  gas  are: 

Conservation  of  Mass:  4^  +  =  0 

at  da 

Conservation  of  Momentum:  _  p 

dt  Sa  ds 

Isentropic  Flow:  -  4^  -  0 


The  last  condition,  with  Initially  uniform  conditions  and  the  equation  of 
state.  Implies  the  following  Isentropic  relations  between  the  gas  proper¬ 
ties: 


(y-0  (y-i) 


The  speed  of  sound  is  found  to  be 


and  can  be  expressed  as 


Rewriting  the  continuity  equation 

Pt  +  P«S  +  U  Pg  >=0 

or  by  introducing  the  acoustic  speed, 

Pf  +  u  Pg  +  p  a’  Og  =  0 

which,  by  adding  and  subtracting  from  the  momentum  equation: 

put  +  p  u  ua  +  Ps  =  0 


yields  two  equations  for  unsteady,  isentropic,  one -dimensional  flow: 
Pt  +  (u  +  a)  Ps  +  pa  [ut  +  (m  +  a)  Us]  =  0 


and 


Pt  +  (u  -  a)Pe  -  pa[ut  +  (tt  -  a)  Usl  =  ® 
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These  relations  are  most  clearly  understood  by  considering  the  "wave 
plane"  or  s t  plane.  By  defining  two  directions  in  this  plane; 


and 


A?' 

dl 


*-'1  =  u  -  a 
It/ 

the  flow  properties  along  these  lines  are  related  by 

d  P  _  —  p  a  du 


and 


dP  =  +  p  a  du 


Integration  of  these  with  the  aid  of  the  isentropic  relations  provides  the 
following: 

and 

(^) 

The  two  directions,  u  +  a  and  a  -  a,  on  the  s-t  plane  are  evidently  the 
paths  of  sound  waves  traveling  respectively  downstream  and  upstream 
through  a  tube  in  which  the  flow  velocity  is  u  .  The  lines  having  slopes 
u  +  a  and  11  -  a  are  termed  downstream  and  upstream  characteristics. 
If  the  velocity  u  and  speed  of  sound  a  are  known  at  a  point  on  a  charac¬ 
teristic,  then  the  relationship  between  u  and  a  is  known  at  all  points 
along  it. 


d  9 


u  +  a  =  const,  corresponds  to  =  u  + 


u  —  a  »  const.  Corresponds  to  =  u  -  a 


APPLICATION  TO  PISTON  MOTION 

It  is  instructive  to  find  the  motion  of  a  frictionless  piston  in  an  in¬ 
finite  tube  of  cross  section  A  having  the  initial  conditions:  (1)  the  piston 
of  mass  m  is  at  rest;  (2)  the  acoustic  speed  behind  the  piston  is  a  ;  and 
(3)  the  tube  is  evacuated  ahead  of  the  piston. 

The  s-t  plane  appears  as ; 
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An  unsteady  flow  region  develops  behind  the  piston  bounded  by  a  sound 
wave  and  the  piston  path.  In  this  region  the  downstream  character- 
isticSj  C  +  ,  connect  the  undisturbed  region  with  the  piston  and  have  the 
property 

y  -  I 

-i— -  u  +  a  -  ac 


It  is  evident  that  the  first  upstream  characteristic,  C„i,is  linear, 
since  along  it  u  =  o  and  a  =  ac  and  ds/dt  =  -ac.  By  considering  the  piston 
motion  to  consist  of  straight  segments  having  discontinuous  changes  in 
slope,  it  can  be  seen  that  an  upstream  characteristic  fan  is  generated  by 
rays  starting  at  each  discontinuity.  The  left-hand  edge  of  this  fan  is  a 
straight  line  by  the  same  reasoning  as  above  --it  bounds  a  region  of  con¬ 
stant  velocity  and  sound  speed.  In  the  limit  of  an  infinitesimal  slope 
change,  the  fan  reduces  to  a  single  characteristic  line  which  is  straight. 
Thus,  the  simple  wave  region  behind  the  piston  is  characterized  by 
straight  diverging  characteristics, C_,  having  the  slope 

=  Up  -  Bp  =  -  Be  +  y  ~  ^  Up 
dt  2 

The  downstream  characteristics,  C^,  on  the  other  hand,  are  curved  gnri 
have  the  slope 

d  »  /  3  -  y  \ 

=  Up  +  ap  =  Be  +  Up 

at  their  intersections  with  the  straight  characteristics. 


The  above  is  the  model  of  an  idealized  gun  having  an  effectively  infinite 
chamber  with  the  same  bore  as  the  barrel.  The  acceleration  of  the  piston 
is 


du_ 

"Tt 


d  3 


PA 


Then,  since  the  characteristics  that  arrive  at  the  rear  of  the  piston  are 
connected  to  a  state  of  known  properties  indicated  by  the  subscript  (c),  the 
following  relation  holds: 


u, 


Bp 


Be 


and  through  the  isentropic  relationship 

2y 


P 


2y  2y 


the  acceleration  becomes 


du. 

“d^ 
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which  can  be,  integrated  to  yield  the  velocity  distance  function: 


P.A 


It  is  usual  to  express  this  in  the  dimensionless  form, 

•  '  »P  ay 


P.A 


9  = 


^0  a  V  2  /  ag 

L  -I 


d, 


which  .can  be  integrated  (Ref.  5)  to  yield 


a  = 


where 


y-l 


P„  A  a 


.  y-l 

U>-^=o 


J±L 

y-l 


y  +  1 


'  a 


and 


SimHarly,  letting 


1  =  t 

m  a_ 


1  == 


dUr 


y  + 1 


[.-ra  =.]’-■  [.-(‘ii);.]’” 


2 

Y  +  1 


These  general  relations  are  shown  in  Fig.  2. 

When  the  launch  tube  is  infinitely  long  or  when  the  projectile  weight 
approaches  zero  (contact  discontinuity),  the  velocity  becomes  ^  _  j  ■ 

The  velocity  of  the  unchambered  gun  having  an  infinitely  long  chamber 
provides  a  convenient  reference  performance  for  comparison  with  more 
realistic  models  to  be  discussed. 

The  effect  of  chambrage  is  discussed  in  detail  in  a  later  section; 
however,  as  a  first  approximation  to  the  effect  of  chambrage,  several 
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authors  (Ref.  4)  consider  the  chamber  to  be  an  infinite  reservoir  from 
which  critical  flow  is  established  immediately  upon  projectile  motion. 
The  energy  equation^ves  for  the  exit  flow; 


•  ^  (^)  - 


and  since 


A, 


O,  e  dj 


U,  = 


J  y  +  1 

from  which  the  downstream  characteristic  value  is  given  by: 


a  +  1^=-^ — I  u  =«  Bj  + 

or  in  a  non-dimensional  form 


a  + 


SO  that 


and 


-  {^) 

=p] 


2y 


y-  1 


(See  Fig.  1) 


Obviously,  this  pressure  can  never  be  greater  than  1.  0  and  the  actual 
pressure  Is  approximately  as  shown  by  the  dashed  curve  in  Fig.  1.  A 
reasonably  good  approximation  to  the  velocity  is  obtained  by  assuming 
that  the  exit  pressure  given  by  the  steady  energy  equation  acts  on  the 
projectile  base  up  to  the  point  where  the  exit  flow  is  sonic. 


SHOCK  COMPRESSION  AHEAD  OF  A  PISTON  AND  FORWARD  FACE  PRESSURE 
PISTON  FORWARD  FACE  PRESSURE 


Free  piston  applications  of  the  gun  principle  utilize  a  compression 
tube  that  is  charged  with  an  initial  pressure  (P,) .  The  forward  face  pres¬ 
sure  on  the  piston  (Pf)  is  derived  in  a  quasi-steady  manner  assuming  it 
is  the  pressure  behind  a  steady  normal  shock.  The  pressure  ratio  across 
the  shock  is  given  in  terms  of  piston  velocity  (Ref.  1): 
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Referring  the  pressures  and  velocities  to  chamber  conditions  Pc  and  ac: 


Figure  3  shows  forward  face  pressure  as  it  varies  with  piston  velocity 
for  a  range  of  acoustic  speed  ratios,  ai . 

The  piston  speed  may  be  computed  from: 

S  =  f  ”  u  da 
•'o  Pr  -  Pf 

for  various  Pi  and  the  appropriate  yc,  /i,  and  ai  . 


SHOCK  COMPRESSION  IN  A  PUMP  TUBE 


When  gas  is  driven  ahead  of  a  relatively  fast  moving  piston,  a  shock 
wave  runs  ahead  of  it.  The  pressure  ratio  across  this  wave  is  related 
to  the  piston  velocity  by: 


yi  (71+  _Z2_  .  ZlZJ: 

Pi  yi +  1 


The  density  and  temperature  ratios  are  given  by: 


^  = 
Pi 


Pi  .  yi  -  1 
Pi  yi +  1 


and 


These  functions  are  shown  in  Fig.  4  for  y  -  1-66 


If  the  end  of  the  piston  tube  is  closed,  the  shock  will  reflect  leaving 
a  state  (3)  of  high  pressure  and  temperature  and  zero  velocity  relative 
to  the  tube.  This  reflected  shock  will  then  reflect  again  from  the  piston 
face  which  la  assumed  to  be  traveling  at  its  original  speed,  resulting  in 
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the  conditions  (4)  (Fig.  5).  If  the  piston  is  infinitely  massive,  shocks 
will  continue  to  be  reflected  from  the  piston  and  end  of  the  compression 
tube.  These  shocks  have  in  common  the  velocity  change  across  the 
shock  equal  to  the  piston  speed,  up.  However,  since  the  sonic  speed 
ahead  of  each  shock  is  continuously  increasing,  the  shock  Mach  number 
becomes  progressively  smaller.  The  shock  compression  approaches  an 
isentropic  process,  and  the  difference  between  shock  and  isentropic  com¬ 
pression  after  the  first  three  shock  transits  is  small.  Figure  6a  shows 
this  in  the  comparison  of  the  pressure -temperature  relation  for  six 
shock  transits  compared  to  three  transits  followed  by  an  isentropic  com¬ 
pression  to  the  same  pressure.  This  leads  to  the  conclusion  that  if  the 
final  pressure  in  this  type  of  compression  tube  can  be  determined 
(empirically  or  otherwise),  the  final  state  of  compressed  gas  will  be 
known.  Practically,  the  final  temperature  may  then  be  related  to  that 
for  isentropic  compression  by: 


Tf 


Tj 


»«11 


y-  1 


T« 

Ti 


JlHi- 

fZi)  y 

V  P,  )  6b) 


This  system  of  heating  of  a  gas  has  application  to  launcher  systems  of 
the  two -stage  free  piston  type. 


EFFECT  OF  CHAMBRAGE 


Chambrage  will  be  defined  as  the  ratio  of  chamber  to  launch  tube 
cross-section  areas.  The  transition  is  assumed  to  be  very  short,  and 
with  increasing  time  the  flow  will  approach  sonic  velocity  at  the  exit. 
Following  Seigal  (Ref.  2)  the  assumptions  below  are  made. 

1.  Flow  in  the  chamber  is  unsteady  isentropic. 

2.  Quasi-steady  conservation  equations  apply  through  the  transition. 

3.  The  characteristics  originate  from  a  rest  state  defined  by  the 
acoustic  velocity,  act ,  at  the  end  of  the  chamber. 

<D].(a  < 

Chamber  Launch  Tube 
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Then,  the  following  equations  apply: 


Unsteady,  one-dimensional  flow  in  chamber 

,  ®ch  “  2  ) 

CpntinuUy  through  the  transition  -- 


(^)  (f:)  ■  (t)  (t)' 


Energy  equation  through  the  transition 


The  above  reduce  to  the  following  relations  between  the  quantities,  u,,  a^. 
11,,  and  a,: 


and 


The  relations  between  characteristic  values,  — )  n  +  a,  and 

characteristic  slopes,  u  ±  a,  across  the  transition  section  Q  and  @ 
can  be  determined.  The  generality  of  the  relationships  is  increased  by 
dividing  all  the  velocities  by  ad,  and  further  normalizing  with  the  Initial 
undisturbed  acoustic  velocity,  a^,  (  bcH  =  ac  until  the  first  upstream  char¬ 
acteristic  arrives  at  the  end  of  the  chamber).  For  use  in  characteristics 
calculations  the  following  functions  are  most  useful: 


Upstream  characteristics  from  chamber  exit  into  chamber 


{^) 


where 


and 
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Downstream  characteristics  from  chamber  exit  into  launch  tube 


The  chamber  exit  velocities 


^  ,  4?-  .  ^  ,  and  ^ 

“ch  “ch  Sch  “ch 


as  functions  of  the  upstream  characteristic  value  from  the  projectile 
path  to  the  chamber  exit. 


“ch 


Graphs  of  these  functions,  (Figs.  7a.  b,  c,  and  d)  permit  deter¬ 
mination  of  the  change  in  the  characteristics  as  they  cross  the  chamber 
exit.  Flow  from  the  chamber  is  limited  by  sonic  conditions  at  the  exit, 
i.  e. ,  u,  =  a, .  Setting  this  limit  and  solving  the  preceeding  equations,  the 
limit  chamber  exit  acoustic  velocity,  is  given  by: 


I 

r 

1  y  + 1 

A.  y  - 1  1 

(  ^  \ 

fi . 

f  2  ^ 

/'■"ch 

iV  1 

1 

1 2(y-i) 

\y  +  i) 

L 

1  y  - 

J 

r 

from  which  the  remaining  velocities,  “jl*  determined. 

The  effect  of  chambrage  on  projectile  velocity  is  illustrated  in  the 
way  the  limit  of  the  downstream  characteristic  value  ~  ^  u,  +  ij  ^ 

varies  with  area  ratio,  •  If  the  chamber  is  infinitely  long,  ach  =  1 ; 
then  the  ultimate  velocity  is  given  by  the  condition  that  ap  =  u 

“  (yrr)  [(^^)  “■L  +  ^i.j 

If  there  is  no  chambrage 

Umax  “  ^  ^  ~  3.03  (y  =  1.66) 

Compared  to 

Umax  ^  1-  1  “c  =  3.33  ac  (y  =  1.66)  (Fig.  7a) 

or  about  10  percent  increase  for  =  4. 

A, 
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When  the  chamber  has  a  finite  length,  ach  is  1.  0  until  the  first  re¬ 
flected  characteristic  reaches  the  exit  after  which  time  it  drops  with  the 
pressure  jit  the  breech  end  of  the  chamber . 

GENERAL  CHARACTERISTICS  SOLUTION  FOR  CHAMBERED  GUNS 

STARTING  POINT- INITIAL  POINT  ON  THE  PROJECTILE  PATH 

In  the  early  stages  of  projectile  motion,  its  velocity  is  small,  and 
the  pressure  acting  on  its  base  may  be  assumed  equal  to  that  at  the  cham¬ 
ber  exit.  In  the  following,  the  launch  tube  ahead  of  the  projectile  is 
assumed  evacuated.  Figure  7c  shows  the  relation  between  a,  and  il,  since 
ac  =  1  at  the  early  time  points.  Graphical  integration  yields  the  coordi¬ 
nates  of  the  first  point: 


Numerical  computations  indicate  that  the  assumption,  a*  =  Sp,  introduces 
negligible  error  for  up  <  0.3.  Beyond  this  point,  characteristics  methods 
of  the  following  sections  are  used. 

Since  «  and  a  are  practically  constant  between  the  projectile  and  exit, 
the  upstream  characteristic  is  a  straight  line  with  the  slope,  up  -  Sp. 

T 


E 


0  3 

At  the  exit  point,  E  : 

Se  =  0 


UE  =  lip  and  ag  =  ap 
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GENERAL  NET  INTERSECTION  POINT 

In  building  up  the  characteristics  network,  new  point,  C  ,  will  be  re 
quired,  which  lies  at  the  intersection  of  a  down  at  re  am  characteristic 
from  point,  A  ,  and  an  upstream  characteristic  from  point,  B . 


^  2  ^ )  "C  +  ac  =  (J 2 

(^)  "C  -  «C  =  (^)  UB  -  Ib 

and  the  local  slopes  are 

uc  +  ac  corresponding  to  the  downstream  characteristic 

“C  -  ac  corresponding  to  the  upstream  characteristic 

The  first  two  equations  are  solved  for  uc  and  ac: 


ac 


"  T  [(^^)  “A  +  aA  -  (^4^)  UB  -  SbJ 
=  [(^)  «A  +  aA  -  acj 


The  coordinates  of  the  point,  C  ,  are  determined  by  extending  the  down¬ 
stream  characteristic  at  the  mean  slope 

-^(uA  +  aA  +  UC  +  ac) 

and  the  upstream  characteristic  at  the  mean  slope 

-^(uB  -  aB  +  uc  -  ac). 

Solving  the  following  equations  simultaneously: 


and 


sc  -  8A  =  (“A  +  ha  +  ITc  +  ac)  dc  -  Ia) 


Sc  “SB  =  -^(uB  -  5b  +  UC  “  ac)  Ctc  -  Tb) 
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The  coordinates  in  the  s  -  1  plane  are  found  to  be: 


»c 


(s  B  -  » A )  +  "tA  (  +  “C  +  °C  ^  -  "tB  ( "B  "  “C  -  ^ 

“A  +  “A  +  2  ac  -  {ug  —  ag) 


and 


®C  “  +  Sa  +  tic  +  ac)  C!c  -  U) 


POINTS  ON  THE  PROJECTILE  PATH 

In  this  case,  a  point,  A ,  on  a  downstream  characteristic  and  a  pre¬ 
ceding  point,  B ,  on  the  projectile  path  will  be  known.  The  next  point,  C  , 
on  the 


projectile  path  is  to  be  determined.  A  method  of  successive  approxi¬ 
mations  must  be  adopted  {Ref.  3).  For  the  general  nth  approximation, 
the  mean  slopes  of  the  characteristic  and  the  projectile  path  are  given  by: 

( [“A  +  aA  +  i«C  (n~l^  +  ac 
=  |8{n)  =  4“  [“B  ^  (n-ll] 

dt^p  2 

The  corresponding  approximate  coordinate  and  velocities  are  given  by: 


(n) 


Afl  ~  9A  +  °  (°)~A  ~  ^ 


o  (n)  —  jS  (n) 


□c  (n)  -  UB  + 


1 


zy  , 

\  r  1 

[ac  ( n  -  1 )]  ^  * 

[ic  (n)  -  fej 

(n) 


=  BA  +  [“A  -  «C  (n)] 
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The  first  approximation  is  made  by  assuming  the  slopes  at  C  are  the 
same  as  at  A  and  B .  That  is 

a  (1)  =  ijA  +  HA 


^  (1)  =  UB 


•Ec  (1) 


”A  +  )  ~A 

“A  +  “A  -  “B 


2y 

«C  (!)  =  UB  +  (ae)^  ^  [fc  (1)  -  "tB] 


ac  (1)=  aA  +  {■^-^)  [“A  -  uc  (1)] 

This  system  converges  satisfactorily  in  three  to  five  steps  and  permits 
calculation  of  the  remaining  coordinate,ac: 

ac  =  Sfl  +  (uB  +  uc)(Tc  -  Tfi ) 


POINTS  AT  THE  BREECH  END  OF  CHAMBER 

The  velocity,  u  ,  vanishes  at  the  wall  resulting  in  the  following: 
ach  =  Downstream  characteristic  value  and  slope,  da/dT 

-  ach  =  Upstream  characteristic  value  and  slope 
A  point,  A  ,  will  be  known  in  an  upstream  characteristic  having  the  value 


then  and  the  point.  B  ,  of  the  breech  end  will  be  given  by: 

-  as  =  “A  -  aA 
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PROCEDURE  FOR  DETERMINING  FLOW  THROUGH  THE  EXIT 

In  general  a  point,  A  ,  on  the  downstream  characteristic  in  the  cham 
ber  and  a  downstream  characteristic,  a,  b,  c,  d,  in  the  launch  tube  are 
known.  This  is  illustrated  below: 


At  the  unknown  point,  C ,  C,  corresponds  to  the  conditions  upstream 
anri  Cj  downstream  at  the  chamber  exit.  Since  the  point,  B  ,  is  not 
known,  it  will  be  necessary  to  use  a  graphical  method.  A  series  of 
values  for  uj  at  C  are  assumed  slightly  greater  than  ua  from  which 
corresponding  times,  Ic »  determined: 


(1) 

(2) 


Tc  =  Ta  - 


2  _ 

"A  +  +  “C  “C 


Determine  the  upstream  characteristic  (^2  )  from 

Figs.  7a  and  d.  —  is  a  function  of 

as  is  .  Multiply  by  Sch  to  obtain 

\  2  /  Bgj,  Bel, 


a,  and  U2  +  a, ;  solve  for  u,  and  a, . 


(3)  On  the  T  -  s  plane,  upstream  characteristics  arriving  at  C 
are  shown  below; 
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The  slopes  of  upstream  characteristics  will  be  (u,  -  a,  +  ir  -  a) 

where  u*and  a  are  determined  at  the  intersection  with  the  down¬ 
stream  characteristic,  a,  b,  c,  d Along  this  characteristic 

~  u  +  a  -=  const.  =  "a  +  aa  j  and  for  the  assumed  up¬ 
stream  characteristics  u  -  a  =  const.  =  * 

u  and  a  at  the  intersections,  a,  ^ ,  and  y,  of  the  characteristics 
may  now  be  solved  for  by  adding  and  subtracting  the  two  simul¬ 
taneous  equations. 

(4)  Plot  upstream  characteristic  value  u  -  a  vs  s 

for  the  intersection  points  above  and  for  the  known  points, 
a,  b,  c,  d,  along  the  known  downstream  characteristic: 


At  the  intersection  read  “»  -  a, ;  from  - —  . 

\  ^  /  \  ^  /  Seh  “eh 

Bp  FiS»  7d>  compute  a, ,  sj  =  ach  —  ^  ^  ^  >  and  tc 

as  was  done  for  the  assumed  values  of  ui . 

(5)  From  Fig.  7a  find  and  compute u,  +  a,. 

5,  and  a,  are  then  determined  at  the  point,  C  ,  on  the  down¬ 
stream  chamber  exit  characteristic. 


Flaw  through  Exit  Prior  to  Arrival  of  First  R»f [acted  Characteristic 
(or  the  Special  Cose  of  the  I  nfinitely  Long  Chamber) 

Until  the  first  reflected  characteristic  reaches  the  chamber  exit, 

the  chamber  downstream  characteristic  value  ii  +  5  =  1  -  Id,.  The 

procedure  is  considerably  simplified  because  the  solution  is  obtained 
directly  from  Fig.  7a  without  successive  trials.  A  point,  B  ,  will  be 
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known  at  which  the  upstream  characteristic,  ub  -  Sb  ,  has  been 

determined. 


The  value  of  “a  -  5^  is  the  same  ub  -  5b  >  and(^^y^)u,  + 

is  found  from  Fig.  7a,  from  which  u*  and  a,  are  computed.  Figure  7d 
gives  the  value  of  ui  and  a,  -  1  -  ^  ^  ^  .  Along  the  upstream  char¬ 

acteristic  (1)-C  in  the  chamber,  u  =•  ui  and  5  =  1,.  The  launch  tube 
characteristic  can  then  be  continued  downstream  from  the  chamber  to  the 
projectile  path  by  previously  discussed  procedures.  The  upstream  char¬ 
acteristic  in  the  chamber  is  straight  until  intersection  with  the  first 
reflected  characteristic. 


PATH  OF  THE  FIRST  REFLECTED  CHARACTERISTIC  IN  CHAMBER 

Throughout  the  chamber  region,  bounded  by  the  first  reflected  char¬ 
acteristic,  the  following  conditions  obtain: 

(1)  The  upstream  characteristics  are  straight  lines  along  which 
u  and  a  are  constant.  The  slope  is  u,  -  1,. 

(2)  The  downstream  characteristics  are  given  by  «  +  a  =  1- 


The  slope  of  the  first  upstream  characteristic  is  -1  since  1  =  0  therefore 

Tch  =  *ch 
1=0 
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The  reflected  characteristic  can  be  completed  as  outlined  for  the  general 
net  intersection  point  with  the  following  simplifications ; 

tic  =  80  ] 


y  at  the  chamber  exit 


u 


c 


UB  J 


a  = 


u 


Computation  Forms 

Figure  8  shows  computation  forms  that  were  found  to  be  convenient 
in  making  the  preceding  calculations.  Figure  8a  was  used  for  Interior 
Net  Points^  8b  for  Projectile  Path  Points^  and  8c  for  Chamber  Exit 
Points.  Using  a  desk  calculator,  50  hours  on  the  average  was  required 
to  complete  a  solution. 


RESULTS  OF  CHARACTERISTICS  CALCULATIONS 
FOR  FINITE  CHAMBER  VOLUME  -  y=  1.66 


Using  the  methods  of  the  preceding  sections,  calculations  were  car¬ 
ried  out  for  helium  (y  =  1.66)  as  a  propellant.  The  following  table  gives 
the  geometrical  variables  considered: 

Chambrage,  A^/A^  Cbambur  Length, 

1  0.5,  1.0,  2.0,  4.0 

4  0.125,  0.25,  0.5,  » 

Any 

The  infinite  chamber  length  case  for  A,/Aj  =  i  is  the  closed  form  solu¬ 
tion  given  in  the  first  section  (Fig.  2). 

Figure  9  is  a  typical  characteristics  net;  A, /A,  =  4  and  xcj,  =  0.5 . 
The  coordinates  are  dimensionless  time  and  distance: 


T  =  ^  8  ^  Pc  Al  a 

■"M  “c  *c 

where 

Pc  ==  Initial  chamber  pressure 
ac  =  Initial  chamber  acoustic  speed 
Al  =  Launch  tube  area 
=  Projectile  mass 
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In  Figs.  10  and  11,  the  dimensionless  velocity  and  time  are  shown  as 
functions  of  dimensionless  distance  for  the  cases  calculated,  u  is  de¬ 
fined  as  uM/ac .  The  reference  solution  for  no  chambrage  and  infinite 
chamber  length  is  included  because  it  provides  a  convenient  basis  for 
normalizing.  It  is  interesting  to  observe  that  the  compensating  effects 
of  chambrage  and  chamber  length  result  in  nearly  the  same  velocity  and 
time  for  equal  chamber  volumes. 


In  Fig.  12  the  above  observation  is  used  to  provide  a  convenient 
summary  of  the  results  of  the  characteristics  calculations.  Here  the 
velocity  is  normalized  to  the  reference  case  (A^/A,  =  1,  xch  -  and 
plotted  against  volume  of  the  chamber  to  volume  of  the  launch  tube  for 
constant  values  of  launch  tube  length,  a  .  It  will  be  noted  that  larger 
chamber  volumes  are  required  for  smaller  ¥ .  This  is  understandable 
when  it  is  considered  that  small  s  is  associated  with  large  acoustic 
speed  and  therefore  rapid  communication  of  the  decaying  chamber  pres¬ 
sure  to  the  projectile  base.  The  apparently  anomalous  results  of  certain 
experimental  launchings  with  cold  and  hot  propellants  are  thereby  ex¬ 
plained.  For  instance,  in  Ref.  6  very  good  correlation  of  launch  velocity 
with  chambrage  was  obtained  using  cold  nitrogen  as  a  propellant.  In 
Ref.  5,  however.  Charters  observes  that  in  high  temperature  helium  the 
velocity  observed  was  best  predicted  by  the  no- chambrage  reference 
formula.  Similar  results  have  been  noted  (Ref.  7)  at  AEDC  for  a  Il^-Ue 
combustion-heated  launcher,  which  has  also  been  operated  with  cold 
helium.  In  these  latter  cases,  the  effects  of  chambrage  and  finite  cham¬ 
ber  volume  approximately  compensate  in  the  hot  firings  ( s  =  4  -  10) ,  and 
in  the  cold  launchings  (i  -  20  -  40)  the  chambrage  affects  the  velocity 

favorably. 


A  criterion  for  selecting  the  volume  of  a  chamber  for  a  given  launch 
tube  is  considered  in  Ref.  2,  which  is  the  arrival  of  the  first  reflected 
characteristic  at  the  end  of  the  launch  tube.  This  leads  to  excessively 
long  chambers  (Fig.  10)  --  from  one-fifth  to  one-half  the  launch  tube 
length.  It  is  suggested  that  a  reasonable  criterion  might  be  derived  fi^m 
Fig.  12  as  the  chamber  volume  required  to  obtain  the  reference  velocity, 
uo  Ui/Aj  -  1,  5rei>  In  this  case  the  chamber  volume  will  be  between  0,1 

and  0.  4  of  the  launch  tube  volume  corresponding  to  i  between  3  and  20, 
which  is  appropriate  for  light-gas  launchers  of  high  performance.  Fig¬ 
ure  13  shows  this  in  simple  form.  It  will  be  noted  that  a  larger  volume 
will  be  required  with  chambrage;  however,  the  chamber  will  be  shorter. 
Such  a  configuration  is  usually  more  convenient  and,  in  addition,  the 
surface-to-volume  ratio  is  more  favorable  with  respect  to  heat  losses. 
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CONCLUSIONS 


A  characteristics  method  is  outlined  for  an  ideal  propellant  by  which 
the  performance  of  a  model  launcher  can  be  computed.  Both  chambrage 
and  finite  chamber  length  are  included.  Dimensionless  time-distance- 
velocity  results  are  presented  for  chambrage  area  ratios  of  1,  4,  and  « 
and  for  several  chamber  lengths  applicable  to  helium  as  a  propellant. 

The  results  are  summarized  to  show  the  effect  of  finite  chamber 
geometry  on  launch  velocity.  The  volume  of  the  chamber  required  for  a 
given  launch  tube  is  presented.  The  influence  of  finite  chamber  volume 
is  particularly  important  in  the  design  of  two-stage  launchers  where  the 
volume  after  compression  may  be  very  small  and  therefore  adversely 
affect  performance. 
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Fig.  3  Pliten  Forward  Face  Protsura  a*  a  Function  of  Speed 
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Temperature  Rotia  vs  Pressure  Ratio  for  Adiabotic  Compressions 
Fig.  6  Adiabatic  Compression  of  Helium 
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c.  Chamber  Exit  Aciwstic  Speed  v«  Flew  Speed 
Fig.  7  Continued 
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d.  Entrance  Velocities,  ul  ond  of 
Fig.  7  Concluded 
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Dimensionlesa  Time 


Olmensionless  Distance,  s 


)«>■ 

CO 


Fig.  9  Typical  Characteristics  Net  {or  Finite  Chamber  Volume 
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Dimensionless  Velocity, 
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Fig.  10  Concluded 
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Chamber  Volume/ Launch  Tube  Volume, 


Fig.  13  ChambM  Volume  Required  to  Produce  Ideal  Launch  Velocity 
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